Tomographic transform on a sphere and topological insulators 
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The tomographic transform was first introduced in the field theory literature long ago. It is 
closely related to Radon transform. In this paper we show how the tomographic transform can be 
implemented on a sphere and apply this result to study surface excitations of a spherical topological 
insulator with a single Dirac cone on the surface. 
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Tomographic transform finds its origin in field the- 
ory in Alan Luther's multidimensional bosonization 
construction^ 1 ^ but explicitly it was written down by 
Charles Sommerfield 3 and he also coined the term. It is 
also closely related to the well known Radon transform 4 
(see also the general reference on this topic §). 

Recently tomographic transform and multidimensional 
bosonization based on this transform" were applied 
in the context of topological insulators^ The general 
reference on topological insulators is (TollTl| . Here we 
will need only the fact that a strong topological insulator 
has an insulating bulk and an odd number massless of 
Dirac fermions on the surface. This can be viewed as a 
defining property of a strong topological insulator. We 
will restrict the discussion to the case of a single Dirac 
cone on the surface. 

The aim of this paper is to implement tomographic 
transform on a sphere and apply the developed formal- 
ism to study surface excitations of a spherical topolog- 
ical insulator. Since the sphere has no boundary, this 
is more appealing, because one does not need to worry 
about boundary conditions which are simply ignored in 
the case when the boundary is infinite plane. It is also 
obvious that condensed matter systems are always of fi- 
nite dimension. 

The relation between quantum spin Hall effect and 
strong topological insulators is well known and discussed 
in detail in the literature. For example, in the topolog- 
ical band theory^— one uses the analogy with quan- 
tum spin Hall effect to count topological invariants of 
3D insulators, and in the topological field theory^ both 
3D topological insulators and quantum spin Hall systems 
are descendants of 4D quantum Hall effect. In this paper 
we show that surface excitations of a spherical topologi- 
cal insulator are the sum of edge states of quantum spin 
Hall disks of the same radius as the sphere. 

The spherical topological insulator was already studied 
in the Refs. [16-18], however our approach is completely 
different and it will turn out that one can still learn some 
interesting facts studying this case. We also briefly dis- 
cuss a strong topological insulator having the topology of 
a torus and make rigorous connection of hydrodynamic 
theory of surface excitations with the topological band 
theory. 

We discuss some technical details first. In the case of 
infinite plane the tomographic transform of a function 



h(x,y) is 

he(0 = J f^V) ' dk J cosk^-e)h{x',y')dx'dy' 

where £ = k ■ x, £' = k ■ x', k = (cos 9, sin 9). It satisfies 
the equation 

J h 2 (x, y)dxdy =JdeJhj(Od( (2) 

where K = {-tt/2 < 9 < vr/2}. This is slightly modified 
version of the tomographic transform. We need to find 
something like this on a sphere. 

The starting point is the analogy with Funk 
transform^ Funk transform of a function / on a 2-sphere 
is defined as 



/(u)ds(u) 



(3) 



where x is a unit vector and the integration is over the 
arclength ds of the great circle C (x) consisting of all vec- 
tors perpendicular to x. This can be viewed as Radon 
transform on a sphere. In the case of the infinite plane 
C (x) are rays perpendicular to two-dimensional vector x. 

This suggests that Euler angles could be the suitable 
choice as parameter space of the tomographic represen- 
tation. Indeed, there are three angles a, /3, 7 needed 
to parametrize rotations in three-dimensions and the 
measure in the space of these angles is given by duj — 
sin (3dadj3d r yi^ Two angles a and j3 can be chosen as the 
direction of x and the third angle 7 as the angle denoting 
the position on the circle C(x). 

Now we pose the problem. Given the function of two 
variables h(9, ip) we need to find a function /i a/ g(7) which 
is real and satisfies the equation 



si 



h 2 {9, V )dVt 



2tt 

sin f3dad(3 I d"fh 2 a p{^) 




(4) 



The pair (a, j3) is an index, 7 is an argument, and u is 
the upper hemisphere. 

To get further insight into the problem we need to 
study the Dirac hamiltonian on a sphere, since multidi- 
mensional bosonization and tomographic transform are 
closely related to each other. According to the Ref. pl| 
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Dirac hamiltonian on a curved surface of a topological 
insulator is 

H = - [V ■ n + n(p x a) + (j> x an)] (5) 

where n is a unit vector normal to the surface and v = 1 
is the electron velocity. On a unit sphere this is (al). 
There is a constant term, |, which was dropped. It will 
be inserted in the final result. 

Then one can observe that the operator (al) and the 
operator a z l z (which is Dirac hamiltonian on a unit cir- 
cle) have the same set of eigenvalues, i.e. integers. 

In (|TJ) we used Fourier analysis. Therefore we can try 
to use spherical harmonic analysis in the case (j4]) and 
expand the function h first 



(6) 
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There are a set of orthogonal functions of three Eu- 
ler angles D l m , m (a, /3, 7) = e lm '^d l m , m ((3)e ma . Their ex- 
plicit form can be found in textbooks,— We will need 
only the orthogonality relation 



D J1 * (a,/3,7)^ 2 , (a, /3,7) 



8^2 



1 



4/1 + 1 

and transformation properties under inversion and com- 
plex conjugation 

d l m , m (7r-P) = (-l) l - m d l _ m , m (l3) (8) 



<> m {fi) - (-i) m - m V m , _ m G8) 



(9) 



These functions are wavefunctions of a symmetric rigid 
rotator. The matrix D 3 gives unitary irreducible repre- 
sentations of rotation group. Therefore one can try to 
use them to construct the tomographic transform. 
Indeed, one can check that the function 



lm 



8tt 2 



+ D l,-m( a + 7T, 7T - /3, 7T - j)hl t . 



(10) 



is real and satisfies Eq.Q. It has the symmetry h a ^{^) — 
p(7T — 7). Further, if fh(8 : <p)dn = 0, i.e. h 
describes surface excitations of some 3D incompressible 
fluid, then J h a p(-~{)d^ = 0, which means that haplj) can 
be viewed as edge excitations of some 2D incompressible 
fluid. 

The next step is to rewrite the Dirac hamiltonian in 
terms of some spinors t^ a p{l)- One can expand the initial 
spinor ip = (ipf ipi) T as 

I in 



Then the hamiltonian becomes 

lm 

Defining new Fermi fields 

ip a pt('y)= a ^n-tD l m , m (a,(3,^) (12) 



I nun' 



VtyM.(7)= a i™iD l m , m+1 (a,f3,"/) (13) 

Imm' 

one can write 

H= j ^(7)(^^(7)dw (14) 



where 
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J+ = J x + iJ y = e ia — -Mcot/3— - i- 

\ap oa smpo7 



9 9 19 

J- = J x -iJ y = e la I -— +icot i 9- 

ap oa sin p 07 



This is due to relations 



J-D l m , m+1 = ^(l-m)(l + m + l)D l m 



J+D l m , m = ^{l-m){l + m + \)D l m , 



m+l 



(15) 



(16) 



J is the space fixed rigid rotator angular momentum 
operator. We want to replace the initial problem of a 
Dirac electron on a sphere with another problem: coher- 
ent states of a symmetric rigid rotator with spin 1/2 (a 
and j3 are polar angles that determine the direction of 
the symmetry axis and 7 describes rotations around this 
axis). This could be done retaining only m! = ±Z, ±(Z — 1) 
in the sums (|12p . (ll3p and multiplying various terms with 
proper coefficients. Now, when one passes to the body 
fixed frame, the action of (crJ) on the rotated spinors 
becomes diagonal (we will do this in the reverse order; 
see (l2"4"|) 'l. To do this one should perform Euler rotation 
?7(a, /3, 7). This is easy to see in the quasiclassical ap- 
proximation when I is large. In this limit the operator 
(<t J) becomes (<rn)Z, because angular momentum is di- 
rected along n, and is diagonalized by rotation to a frame 
in which z axis points towards n. Under rotations of the 
coordinate frame, the spin part of the wave function and 
its angular part transform independently. The inverse of 
the operator which acts on spin part is 



\ sin § • e ^ a ~^/ 2 



-i(a-7)/2 N 



cos f • e l{a+ '<^ 2 



(17) 
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The transformation of angular part in the sector with 
momentum I is described by the D-function 

ll>lm = Dln'm^lm' ( 18 ) 
m' 

There is the following formula 

U^ajj = (na) (19) 

where n = (sin /3 cos a, sin /? sin a, cos /3), which will be 
used later. 

Spinors in the rotated frame i/i^pi"/) can be bosonized 
introducing boson fields 

€h) = v^E e -w i bi a yy h + h - c A ( 20 ) 

l>0 v ' 



compute the derivative id/dt. Due to simple equations 
of motions for cf>^ and <p^ 

d t <j) {1) + <9 7 </> (1) = 0, d t (j) {2) - <9 7 (2) = (25) 

one can replace this derivative by ^fid/d-y and then inte- 
grate by parts, d/d-y now is multiplication by a constant. 
On the other hand, action of (at) on rhs of can be 
replaced by — (a J), J acting on D-functions only, not 
on U or bosonic exponents (this is due to formulas ([15)) 
and |H])). Using (aJ)fj- 1 = JJ-^aP) (J does not act 
on U) where P is the body fixed rigid rotator angular 
momentum operator 

/ d d 1 d \ 



6 (2) 



( 7 ) = V2^^ 7r [6g t (0e^ + /i.c] 
i>o 



(21) 



e is needed for regularization of the theory and one should 
take the limit e — > in the final result. Boson creation 
and annihilation operators satisfy the commutation rela- 
tions 



[b ( * ) (l),b^(j)]=S ik Si j S nta , 



(22) 



and their time dependence is e . Correlators of bosonic 
exponents are 

(exp[-0( 1 )(i, 7 )]exp[^ 1 , ) (O, 7 ')]> 



h _ e -i(t-7+70- e ] 
oo 



(23) 



1=0 



and for one should change the sign of 7 and 7'. 
The usual Fermi field is given by 



2^ 



X < e -('-7)/2 e (1) 



+ e -«(*+7)/2 e (2) 



-i,-r. 




_J_ it 



V2TT1 V A-i,- 



D l,- m J V2TTT 



2 + 1. -m 

21D 1 



l,—m 

(24) 



Here i7 stands for U(a, /3, 7), the arguments of the D- 
functions are a, /3, 7, i s the shortening for 0^0(7) 
and di7 = sin /3dadf3. This expression seems complicated, 
but it has a clear meaning. We will demonstrate this 
by deriving equations of motion for this field. First, we 



P-=Pi 



d n d Id" 

— -zcot/?— — 
op 07 sin p oa / 



07 



which has anomalous commutation relations [Pi , Pj ] 
—ieijkPk and acts on the D-functions according to 



p n l 

r — 1J m'-\-1,m 



^(l-m')(l + m' + l)D l m , 



P+ D rn'm 



y/ (I - m')(l + m' + 1)D 



m +l,m 



P D 



and the fact that spinors in square brackets are eigen- 
functions of (<tP) with eigenvalues —I and I + 1, one can 
show that 



(— — 



This is correct equation of motion. 

Now we will outline the derivation pairwise correlation 
functions (ipip ) and (tp 'ijj) from (|24l) (details of the cal- 
culations can be found in appendix [AJ . Using (1231) one 
can see that integration over 7 and 7' in the fermionic 
correlator picks up the sector with momentum I from the 
bosonic correlator. First of the two spinors in the square 
brackets in (|2"4"]l are important for the calculation of the 
correlator (ip4>' t ), and the second ones for the correla- 
tor (ip^ip) only. There are two independent contributions 
to the fermionic correlator emerging from two different 
bosonic exponents. The part obtained from the terms 
in the second quare brackets in (|2"4l can be converted 
to the the part obtained from the first square brackets, 
but where the integration is over the lower hemisphere, 
by using transformation properties of D-functions 



4 



Calculating the emerging integrals of the product of three 
D-functions (cos/3 = <1q O (0), sin/3 = \/2d\ (P)) which 
are now over the entire sphere, using formulas in [20l ]. 
it can be shown that one obtains correct expressions for 
fermionic correlators. For example 

(^ t (M,*#j(MV)) 

= £ e -«(m/ a )£+il+i yjm(ff) m * mi e>, v >) 

lm 

which coincides with the expression obtained by conven- 
tional means using eigenfunctions of the operator (al) 
with eigenvalue I 

1 1 + 771+1 /l\ ll-m /()\ , . 

V^Ti-W^ + V^TT(iJ y '- +1 (26) 

and eigenvalue —I — 1 

ll-m fl\ T _ ll + m + 1 /0\ 

The necessity of introducing additional factors 
+ 1) in Ea. (l24l) can be seen by explicit cal- 
culations (see appendix [XJ. 

The operators do not anticommute. Therefore one 
should insert Klein factors in (|24D that ensure anticom- 
mutation relations. In the continuum casei one uses the 
limit k — > of some boson creation operators for this 
purpose, where k is momentum. This is not possible in 
the discrete case we are dealing with. We will introduce 
some finite but very small compressibility of the liquid 
and static bosonic field a instead. Then if one adds a 
term b^g/S 1 ^ 2 to and inserts the operator 

Q n = e ™£ n -< n [# n ) +&& , ]« 1/2 

into the integrand in (j24|) . then ((24)) will satisfy the re- 
quired anticommutation relations 

{^j(t ) x),V'k(*,x')} = SikS(x - x) 

We assume that 6 is small. Introducing the operator b^L 
into the theory modifies the electronic density operator: 
there will be an extra term proportional to 5 1 / 2 which 
tends to in the limit 8 —¥ 0. This is consistent with the 
incompressibility of the liquid. 

One can trace the analogy between the scheme devel- 
oped in this paper and that of Luther's almost in every 
step. Luther also showed that bosonic exponents in his 
scheme lead to correct expressions for general correlation 
functions. His analysis relies only on the two-point corre- 
lation functions, equations of motion for the fermi fields 
and anticommutation relations. One can directly extend 
this analysis to the present case and there is no need to 
repeat them here. 

The formula (|2"4"]l has no useful applications. It was 
derived only to show that bosonization on a sphere can 
be consistently carried out. 



Though we mentioned the relation of our analysis 
to topological insulators throughout the paper several 
times, it was quite abstract. Now we will make this 
relation more concrete. Following the Ref.[8( one can 
describe low lying excitations of a strong topological in- 
sulator as surface deformations of a two-component 3D 
incompressible liquid confined by a smooth potential well 
with the Hamiltonian 

H = ^ Pa eER 2 J [h\{6, if) + h 2 2 (9, <p)] dSl (28) 

where E is the electric field of the confining potential on 
the surface, po is the density of the 3D electronic liquid, 
J hi(8,(p)dCl = 0, i = 1,2 and R is the radius of the 
sphere. This theory is a modification of the hydrody- 
namic theory of the edge excitations of fractional quan- 
tum Hall systems. 22 Using (HJ) and (jTUJ) one can rewrite 
(|2"5f as follows 

H=±neERj dilJ [h 2 laP (-y) + hj aB ( 7 )] d 1 (29) 

where hi a p can be viewed as surface deformations of a 
2D incompressible fluid with the density n — Cp^R and 
having a shape of a disk of radius R, C is a numeric con- 
stant {h a p has been rescaled by a factor C -1 / 2 compared 
to (fTOl) '). We have assumed in (|29| that electric field of 
the potential which confines this liquid is E. In analogy 
with the paper Q we assume that hi a p (7) are edge states 
of quantum spin Hall system, i.e. they have equations of 
motion 

d t hi a (i(j) = {-lyvd^hiap^) (30) 

where v — eE/hn. Thus the constant C is not arbitrary: 
C = eE/hp vR. The theory ([H can be easily 

quantized 

H= f dfi£ E wi6 $ + (0 6 $(0 (31) 

Z>0 i=l,2 

where the boson operators b^L were introduced earlier 

mim. 

Thus we have found that surface excitations of a spher- 
ical topological insulator with a single Dirac cone on the 
surface are the sum of quantum spin Hall edge states. 
Since the transform (|10[) is non-local, every edge state 
that has a fixed location in the tomographic representa- 
tion is spread over the entire surface of the topological 
insulator. 

Now we will briefly discuss a toric topological insu- 
lator for which one can make a direct connection with 
the topological band theory. This case is simpler than a 
spherical case. 

Suppose that there is an ideal crystal of a strong topo- 
logical insulator with atomic spacing a having » 1 
unit cells along x and y directions and periodic boundary 
conditions imposed in these directions, z = is one of the 
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surfaces of this crystal. The function h(x, y) defined on 
this surface and having the property J h(x,y)dxdy = 
can be expanded as (L = Na) 

h(x, y) = J2 h mn e 2 ™ {mx+ny)/L (32) 

n in 

and there can be denned a new function 

oo 
2 = 1 

+h- lm ,- ln e- 2 ^ L ^] (33) 

where 

a 

and (mn) denotes a pair of coprime integers m and n 
such that to > and n arbitrary or one of the pairs 
?7i = 0,77 = 1 or 777 = 1,77 = 0. (77777) are in fact Miller 
indices and / has the meaning of a winding number. It is 
easy to see now that 

J h 2 (x,y)dxdy = L^ f imn> h 2 (mn) {Qdt (34) 

(mn) 

where the sum is over all possible pairs (77777) . This means 
that we have rewritten the 2D hamiltonian as a sum of 
ID hamiltonians over closed loops of length Lr mn y 

One can see that if h(x,y) describes surface excita- 
tions of 3D incompressible fluid, h o = 0, then /i( TOn )(£) 
describes edge excitations of a 2D incompressible fluid, 

The above procedure applies also to 2D fermionic 
hamiltonian 

m.ln 

(35) 

l,{mn) 

The operator ma x + na y can be diagonalized by a ro- 
tation. Diagonalized hamiltonian is 2nb z v\J m 2 + n 2 /L. 
Eq. (j35|) describes surface excitations of a strong 
topological insulator for small momenta ly m 2 + n 2 <C 
N. Surface deformations of two incompressible liquids 
/ij( mn )(t, £), i = 1,2 will satisfy equations of motion 

<9t/i,( TO „)(i,0 = (-l) l i%/7 i ( m „)(i,£) 

These equations are consistent with the spectrum 
2-7T ■vlyrn 2 + n 2 / L of the fermionic hamiltonian (I35[) and 
the exponent e ±27r2 ^/ i ' m ™) we have chosen in (|3"3"|) . It 
is clear that one can construct bosonization scheme in 
analogy with the Ref. [l|. 



Now it is easy to make parallel with the Z2 invari- 
ants of 2D band structures^ Let H(k x ,k y ,k z ) be a 
tight binding Hamiltonian of a strong topological insu- 
lator (see, e.g. the Ref. [Ill) and Hr mn )(2Trl/Lf mn j, k z ) = 
H(2irlm/L, 2irln/L, k z ). There is one to one cor- 
respondence between H ^ mn ^{2-Kl / L ^ mn y z) and the 
edge state hamiltonians (|35D for fixed direction 
(mn). The hamiltonian if( mn )(27r/yL( m „), k z ) is pe- 
riodic with the period ^ - along k^ direction. 
This means that the 2D lattice has Brillouin zone 
( 5— < ke < —Z—.-Z < fc z < il\ it follows from 

the topological band theory that all hamiltonians H(jn n ) 
for fixed (77777) are hamiltonians supporting quantum spin 
Hall effect. The easiest way to see this is to use the sim- 
ple countin g argu ment due to Roy2^, or the approach of 
the papers |12ll25l |. In other words Zi invariants of the 
planes (77777) are all equal to 1 in the case of a strong 
topological insulator. This means that the fields h i ^ mn ^ 
i = 1, 2 (in the hydrodynamic theory) are the edge states 
of quantum spin Hall systems. 

There were some difficulties in the continuous casei* 
It is clear now that these difficulties can be overcome by 
considering a discrete version of the tomographic trans- 
form. 

The assumption of equal number of sites along x and 
y directions is not essential. If there are different num- 
ber of sites, then this would only lead to unnecessary 
complication of the analysis. 

In summary, we have shown how the tomographic 
transform can be implemented on a sphere. Using 
this transform, we constructed Luther's version of the 
bosonization on a sphere and showed that surface exci- 
tations of a spherical topological insulator with a sin- 
gle Dirac cone on the surface are the sum of quantum 
spin Hall edge states. However, this is not true for the 
entire topological insulator: strong topological insulator 
can not be presented as a sum of two-dimensional topo- 
logical insulators. We also constructed discrete version 
of the tomographic transform on a torus and rigorously 
showed how the connection can be made between hydro- 
dynamic theory of surface excitations of a non-interacting 
topological insulator and the topological band theory. 
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Appendix A: Calculation of correlation functions 



In this appendix, we will show how to compute two-point correlation functions starting from the Eq. (fMf . After 
calculation of bosonic exponents and integration over 7 and 7' one obtains 



lm\m2 



d ul±p^D l _ u _ mi (a, (3, 0)D l * h _ m2 (a, (3, 0) + j dd-^D{_ mx (a, /3, 0)D\*_ m2 (a, (3, 0) 
Using Eqs.© and © it can be shown that the second integral in square brackets equals 



±^D l _ h _ mi (a, (3, 0)D l * lt _ m2 (a, f3, 0) 



where l U means that integration is over the lower hemisphere. The integral 

dn l±pL D l_ l _ mi (a> j9j 0)D l * lt _ m2 (a, (3, 0) 
can be calculated using the following integral of the product of three D-functions 



D J1 , 



duj 



, (uj)D 32 , HL> 33 , M — 

\rn-i V / m' 2 m 2 y > m' 3 m 3 \ ' g^2 



Ji 32 33 \ 3i 32 33 



m'i m' 2 m' 3 J \mi ni2 1TI3 



(Al) 



and the fact that 



l±COS/3 x 1 . ! 

2 ^^iiC/ 3 ): -^sin/3 = d 10 (^) 



/ 71 79 73 \ 

Here I ^ ^ ^ J is a 3j-symbol which is non-zero only when 7713 = —mi — 7712. For ji = 22 an d 33 = 1, 3j'-symbols 
can be easily calculated using tables ir*2£. For example, 



/ I 1\ = t _ m 2m / J J I 

m -m 0y 1 J [2Z(2Z + 1)(2Z + 2)] 1 / 2 ' -m - 1 



= (-!)« 



2(Z-m)(Z + m + l) 
2Z(2Z + 1)(2Z + 2) 



1/2 



The result coincides with the one obtained by using the eigenfunction (J2HJ) of the operator (al). 
For the other correlator (ip^ip), after integration over 7 and 7', one has 



x / dn 



(ft(0,()>, ( p')i>(t,9, V )) = ±- £ ^^^(-l)" 11 -" 12 ^^^)^^^') (A2) 
2JI>{ i _ mi (a, 0, 0)£^_ m2 (a, /3, 0) + D\_ h _ mx (a, /?, Q)D{* 1; _ TO2 (a, /3, 0)i^ 



(A3) 



+a+(a, (3W2lDl_ mi (a, (3, 0)Dj* 1; _ m2 («, ft 0) + (a, /3) V2/£>|_ 1 ,_ mi (a, (3, 0)D\*_ m2 (a, /?, 0) 
DU+i,- mi («, ft 0)Z?i! /+ i,- roa («> A °)^^ + («, ft 0)^,_ m2 (a, /3, 0)i^ 

+a + (a, P)V2W l _ l+lt _ mi (a, 0, $)D-l,- m . («, ft 0) + <r_ (a, (3)V2lD l _ h _ mi (a, f3, Q)D l * l+li _ m2 (a, (3, 0) 
where 



or explicitly 



. , « ~t I m 1 f sin/3 e~ M (l + cos/3) 
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Let A be the expression in the first square brackets in (| A2|) and B in the second. Then using the Eqs. © and (0) it 
can be shown that 



/ Adn = [ Bdii 

JL Ju 



The integral J AdVl can be calculated using (|A1|) . For example, 

AndQ = — -Sr, 



and 

<$(0, 9\ tp'^t, 0, <p)) = ]T ^e^ +1 ^Y lm (0, v)Y? m (0',<p') 

lm 

which coincides with the expression obtained from (|27[) by conventional means. 
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